Fluctuation Probes of Early-Time Correlations in Nuclear Collisions 
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Correlation measurements imply that anisotropic flow in nuclear collisions includes a novel tri- 
angular component along with the more familiar elliptic-flow contribution. Triangular flow has 
been attributed to event-wise fluctuations in the initial shape of the collision volume. We ask two 
questions: 1) How do these shape fluctuations impact other event-by-event observables? 2) Can 
we disentangle fundamental information on the early time fluctuations from the complex flow that 
results? We study correlation and fluctuation observables in a framework in which flux tubes in 
an early Glasma stage later produce hydrodynamic flow. Calculated multiplicity and transverse 
momentum fluctuations are in excellent agreement with data from 62.4 GeV Au+Au up to 2.76 
TeV Pb+Pb. 

PACS numbers: 25.75.Gz, 25.75.Ld, 12.38.Mh, 25.75.-q 



I. INTRODUCTION 

Measurements of two particle correlations in nuclear 
collisions exhibit a complex pattern of ridges, bumps, 
and valleys as functions of relative pseudorapidity A77 = 
Vi ~ V2 an d azimuthal angle A</> = 4>i — (f>2 [IHS] ■ I n & re- 
cent advance, Alver and Roland showed that much of the 
A(j) dependence of these correlations can be described in 
terms of anisotropic flow |10j . Their key realization was 
that a novel triangular flow contribution ~ W3 cos 3A0 
is needed to describe the data. Such contributions have 
been attributed to fluctuations of the geometric shape 
of the collision volume from event to event [TDT - BTj . A 
unique shape is determined in the first instants of each 
collision event as the nuclei crash through one another. 

Interestingly, fluctuations have been of measured in 
nuclear collisions for many years, but for different rea- 
sons. Event-by-event fluctuations of the multiplicity, 
mean transverse momentum, and net charge are viewed 
as probes of the QCD phase transition, although signa- 
ture behavior has yet to be seen [221424) . Such observa- 
tions are explained by the variation of density, geometry, 
thermalization, and flow; see, e.g., [25H29] . 

In this paper we explore the common influence that 
the early-time dynamics of the system has on correla- 
tions, flow, and fluctuations. We begin by asking how 
geometric fluctuations impact other event-by-event ob- 
servables. In Sec. |TlJ we recall the relationship between 
correlations and anisotropic flow. Next, we discuss the 
relationship between fluctuations and correlations in Sec. 
Ill We exploit integral relationships obtained in Ref. [3U] 
to marry the dynamical language of flow and correlations 
to the statistical formulation of fluctuation observables. 
This yields model independent properties of experimental 
observables. In particular, we find that geometric shape 
fluctuations alone cannot explain measured multiplicity 
and transverse momentum fluctuations. 

We next ask what fluctuation and correlation measure- 
ments can reveal about the early-time dynamics of the 



collision system. The first evidence that correlations orig- 
inate at early times in the collision is the long rapidity 
range of the ridge [3TJ [32] ■ Correlations show a ridge-like 
peak near A(f> = 0. This ridge extends over a broad range 
in A?7, as do away-side features centered near A</> ~ it [5]. 
Causality dictates that correlations over several rapidity 
units must originate at the earliest stages of the collision 
[5T| 152] . With that in mind, we described the ridge as 
a consequence of particle production in an early Glasma 
stage followed by transverse flow in Refs. [5TH54"] . Our 
description is part of a broader family of models in which 
particles are initially correlated at the point of produc- 
tion mi in GEM- 

Identifying the impact of anisotropic flow on these cor- 
relations adds considerable credence to this observation 
[Tl] . Anisotropic flow is well understood as an early time 
effect, since it is generated in part by the geometric con- 
figuration of participant nucleons in the colliding nuclei; 
see, e.g., [SI 02]. Correspondingly, the measured v 2 co- 
efficients vary little with rapidity. 

To illustrate how correlations, flow, and fluctuations 
result from the early-time dynamics, we apply the gen- 
eral framework of Refs. [321 [33] in Sec. |IV| We obtain 
expressions for the correlation function and its Fourier 
coefficients, |23l, (24 1, and (25). This provides a unify- 



ing framework for understanding hydrodynamic v n stud- 
ies together with earlier work on the ridge. We also derive 
expressions for transverse momentum fluctuations, (26) 



and ([27)). 

In Sec. [VJ we use fluctuation measurements to extract 
information on the particle production mechanism. We 
focus on the Color Glass Condensate formulation in Ref. 
[31H331 43 , and argue that dynamic multiplicity fluctu- 
ations 1Z can provide information on particle production 
that is independent of flow. Transverse momentum fluc- 
tuations (SptiSpt2) provide similar information, although 
the results are somewhat sensitive to radial flow (but not 
the v n ). Constraining the flow contribution by calculat- 
ing (p t ) and V2, we compute the fluctuations for 1Z and 



(8pti8pt2). We find that the same model that described 
the energy, target-mass, and p t dependence of the ridge 
also describes transverse momentum fluctuations mea- 
sured at the Brookhaven Relativistic Heavy Ion Collider, 
RHIC, and the CERN Large Hadron Collider LHC. 



II. 



CORRELATIONS AND FLOW 



Correlation measurements commonly center on the 
pair distribution 



MPl,P2) 



dN 



dr]id 2 ptidrj 2 d 2 pt2' 



(1) 



where r\i is pseudorapidity and pu — \pu\ is transverse 
momentum of particles i = 1,2. In the absence of cor- 
relations, p2(pi,P2) -> Pi(Pi)pi(P2), where the single 
particle distribution is pi(p) = dN/di]d 2 p t . Experiments 
typically report ratios of p 2 to either p\p\ or mixed event 
pairs integrated over ranges of p t . To keep the notation 
simple, we will not make the r\ dependence explicit unless 
needed. 

To illustrate the connection between flow and correla- 
tions, recall for the moment the traditional picture of flow 
in which event-by-event fluctuations are neglected. It has 
long been known that collisions at non-zero impact pa- 
rameter b produce anisotropic flow |41[ I42j . Anisotropy 
derives from the change in the shape of the collision vol- 
ume with respect to the reaction plane, i.e., the plane 
spanned by b and the beam direction. The distribution 
with respect to this plane is 



Pi{VfA RP ) =Pi{l' 



oo 



2 ^ v n (p t ) cos[n(4> - 



•V>«p)]}, (2) 



where the coefficients depend only on the magnitude 
Pt = \pt\ and r\. An experimental analysis that does 
not identify the reaction plane measures the distribution 
(27r) _1 J pi(p t ;ip RP )dip RP =Px, where the "bar" denotes 
average over </>. The second equality follows because p\ 
is a function of — ip RP . 

Anisotropic flow introduces correlations because pairs 
from the same collision event have the same reaction 
plane. The reaction plane averaged pair distribution can 
only be a function of the relative angle A0, so that 



OG 



P2 = P 2 {! + 2 / _, a n {Pi,P2) cos(nA</>)}. 



(3) 



If geometry is the only source of correlations then 
P2{Pti,Pt2) is the product p 1 (p tl ;ip RP )p 1 (p t2 ; ip RP ) av- 
eraged over ip RP . In that case the Fourier coefficients of 
(pi) and (pi) are related. For n = we find 



P 2 {Pl,P2) -* Pl(Pl)Pl(j>2)- 

while for n > 2 



uncorrelated (4) 



a n (pi,P2) -> v n (pi)v n (p 2 )] uncorrelated (5) 



these results hold only when fluctuations may be ne- 
glected. Momentum conservation contributes to oi, mod- 
ifying the Vi(pi)vi(p 2 ) term, as discussed by Borghini 
and others [ifHIg] . 

Fluctuations introduce further anisotropy because the 
shape of the collision volume is different in each collision 
event. In collisions of identical nuclei, the event-averaged 
interaction volume is symmetric in <fi and fixed by b and 
^ RP . If all events of a given b had the same interaction 
volume, we then would expect only even n to contribute 
to ([2]). Shape fluctuations cause the flow parameters v n 
to vary from event to event and allow odd n to contribute. 
The average pair distribution has the same form as (pL 
but with coefficients 



a n {Pi,P2) = (v n (Pi)v„(p 2 )} , 



(6) 



where the brackets denote average over events (includ- 
ing all event shapes) [T^]. The measured azimuthal de- 
pendence of two particle correlations are reasonably de- 
scribed by (pi with n = 1,2 and 3. We emphasize that 
fluctuations in shape alone cannot alter p 2 - this requires 
further dynamical fluctuations which we discuss in the 
next sections. 

We remark in passing that one often discusses the 
shape fluctuations in terms of a geometric eccentric- 
ity e„. If one assumes that the relation between e„ 
and the resulting anisotropy of the fluid flow v n is ap- 
proximately deterministic, then fluctuations of the ratio 
v n {pi)/e n are negligible. In that case (v n (pi)v n (p 2 )) w 
{v n (pi)/£n)(vn{p2)/e n )(£n)- This factorization conjec- 
ture seems plausible, but currently requires further the- 
oretical investigation |4"9] . 



III. 



FLUCTUATIONS 



Fluctuation measurements study the variation of bulk 
quantities, such as multiplicity or average transverse mo- 
mentum, over an ensemble of events [22, 50, 51 . Suppose 
that fluctuations of multiplicity N result in a variance 
a N = (^ 2 ) — (^0 2 - Uncorrelated particles would be de- 
scribed by Poisson statistics, for which a 2 — > a 2 — 
(N) . Correlations give rise to a difference a 2 , — o 
which we characterize by the dynamic variance 



stat 
2 

stat ' 



n 



(N 2 ) - (N) 2 - (N) 



(7) 



as discussed in [3U]. Similarly, many describe the dy- 
namic fluctuations of transverse momentum using the 
covariance 



(SptiSp t2 ) 



T,i*4 SptiSp-. 



w 



(N(N-l)) 



(8) 



where Spu = pu — (pt) and the average transverse mo- 
mentum is (p t ) = (P t )/(N) for P t = YuiPti the to- 
tal momentum in an event [52H55J . This quantity van- 
ishes when particles i and j are uncorrelated. Note 



that one can write (5ptiSpt2) in terms of the variance 
(Xp = ((P t — N(p t )) 2 ). In the absence of correlations, 



IV. 



SOURCE OF FLUCTUATIONS 



that variance is up stat 



= <i\o((p?) - (p t ) 2 ). o 



ne can 



show that (SptiSpt2) is the difference ap t — <J PtStat di- 
vided by the average number of pairs (N(N — 1)). 

To see the connection between fluctuation and corre- 
lation measurements, observe that the integral of the 
pair distribution p 2 gives the average number of pairs 
(N(N - 1)). We then write Q as 

1 



K = 



where we define 



{NY 



-(pi,p 2 )dpidp 2 , 



(9) 



KPl,P2) =P2(Pl,P2)-/Ol(pi)pi(p 2 )- (10) 

Similarly, we can write pi in terms of the correlation 



function ( 10 ) as 



(Sp t iSp t2 ) 



«Pl«P2 /Ar , Ar ^5ptl5p t 2- (11) 



(N(N-l)) 

We stress that the densities p\ and p 2 are event-averaged 
quantities. 

Fluctuation measurements probe the overall strength 
of correlations in a manner that is independent of the 
anisotropic flow. To see this, we combine fcfy with ( 10 ) to 
find that the cos nA<p contributions vanish on integration 
over < A(f) < 2tt. We obtain 



n = 



i 



(Ny 



"(pti,Pt2)dpidp 2 , 



which depends only on the (j> averaged function 
r{pn,Pn) = P 2 (pti,Pt2) - Pi(pti)~Pi(pt2)- 



(12) 



(13) 



It is easy to understand why 1Z is not sensitive anisotropy 
- N simply counts particles irrespective of where they are 
flowing. Similarly the 



(Sp tl Sp t 2) = 



, , r(p tl ,p t2 ) , . 

dpidp2 (N(N~l)) SPtlSPt2 > (M) 



Here too the cos nA<p contributions to fcfy vanish on inte- 
gration. These fluctuations are independent of <fi because 
our definition of p t i disregards direction. 

Equations ( 12 1 and ( 13 ) have two striking implications 



when combined with flow results from the previous sec- 
tion. First, if the variation of the initial geometric shape 
of the collision volume is the only source of fluctuations 
then Q implies that 1Z and (dpndptz) must both van- 
ish. Experiments have measured these quantities and 
they are both nonzero, as shown in Sec. |Vj Second, the 
amplitudes r and p 2 determine the overall magnitude of 
correlations, as we see from ([3]). Anisotropic flow and mo- 
mentum conservation determine the coefficients (pi) and, 
therefore, the relative height of the near-side "ridge" at 
Acf> = and away-side features near A</> = it. However, 
interpretation of the evolution of the ridge height with 
beam energy or centrality requires an understanding of r 
or TZ. 



Nuclear collisions vary sharply from event to event due 
to differences in the number and configuration of the nu- 
cleons struck in the initial impact. Each strike adds to a 
transient color field that lasts a proper time of roughly 
To ~ 1 fm. This field comprises an array of flux tubes con- 
necting the fragments of the highly Lorentz-contracted 
nuclei along the beam direction. The number of partic- 
ipants determines the color charge and thus the overall 
strength of the fields. The flux tubes fragment after To, 
driving soft particle production. We emphasize that flux 
tubes arise naturally in QCD and have long been the 
core of phenomenological models such as PYTHIA. In 
the next section we will focus on the Color Glass Con- 
densate description, which incorporates these features in 
the high density environment produced by nuclear colli- 
sions and allows for systematic computations. For now, 
we keep the discussion more general. 

We assume that the number and geometrical distribu- 
tion of flux tubes is the most important source of fluctu- 
ations. Their fragmentation leads to density correlations 
described by 



c(xi,x 2 ) = ti 2 (xi,x 2 ) -ni(xi)m(x 2 ), 



(15) 



where n\ and n 2 are the single and pair densities. In 
the absence of correlations, n 2 (xi,x 2 ) — > ni(xi)ni(x 2 ) 
so that c vanishes. The integral of n 2 over position gives 
the number of pairs averaged over events, (N(N — Vj), 
so that the integral of c is 1Z(N) 2 , with 1Z given by (|7|). 
We now take the flux tubes to be longitudinally boost 
invariant, so that the correlation function only depends 
on transverse coordinates r t = r lt — r 2t as well as the 
average R t = (ri^+r 2it )/2. We further assume particles 
from the same flux tube are primordially correlated and 
that correlations with the reaction plane arise due to the 
distribution of tubes. Since the transverse size of the flux 
tube is small, the primordial correlations reflect common 
spatial origins. We then write 



c( Xl ,x 2 )-^(iV) 2 J(r t )p Fr (R t 



(16) 



Here, p FT (Rt) is the probability distribution for finding 
a flux tube at a transverse position R t in the collision 
volume. This function describes the distribution of shape 
fluctuations discussed earlier. Following Ref. [22 [S3], we 
take p FT to roughly follow the participant distribution of 
the colliding nuclei 



p FT (R t ) 



nR\ 



1 



3. 



(17) 



for R t < Ra : and zero otherwise. 

We now discuss the impact of these long range correla- 
tions on the final-state distributions. Comoving partons 
locally thermalize as the flux tubes fragment. Pressure 
builds and transverse expansion begins. In this process, 
partons from a flux tube at an initial space-time point 



(i, x) will eventually acquire a final flow four velocity u^. 
In a blast wave model u^ is taken to have a Hubble-like 
correlation, while hydrodynamic calculations provide a 
more realistic u^. In either case, the Cooper- Frye single 
particle distribution is 



Pi(p)= //(x,p)dT, 



(18) 



where / is one-body phase space distribution function 
and dT = p^da^ is the element of flux through the 
four dimensional freeze out surface. In local equilib- 
rium the distribution has the Boltzmann form /(x, p) = 
(27r)~ 3 exp{— pVup/T}. The temperature T and fluid 
four-velocity u^ are generally fixed by hydrodynamics, 
which enforces the local conservation laws. In keep- 
ing with the boost-invariant distribution (16), we as- 



sume that freeze out occurs at a proper time tf, so that 
p^da^ = TFirit cosh(y ~rj)drjd 2 rt, where r\ = (1/2) ln((t+ 
z)/(t — z)) is the spatial rapidity. 

The pair distribution has an analogous Cooper-Frye 
form 



Pa(Pi,P2)= / / 2 (xi,pi,x 2 ,p 2 )dri(ir2, 



(19) 



where f 2 is the two-particle Boltzmann distribution func- 
tion. In local thermal equilibrium the two particle distri- 
bution is 



h = ™ 2 (xi,x 2 ) 



/(xi,Pi)/(x 2 ,p 2 ) 
n(xi) n(x 2 ) 



(20) 



where n\ and n 2 are the single-particle and pair densities 
discussed earlier, with n(x) — J /(x, p)d 3 p. This expres- 
sion satisfies a generalization of the Boltzmann transport 
equation for / 2 ; the factors of / cause the generalized 
collision terms to vanish just as they do in the one body 
equation. We omit momentum and energy conservation 
terms that do not contribute to 1Z and (Sp t iSp t2 ). 

We understand (20) as follows. In local equilibrium 



we can divide the system into fluid cells, each of which is 
in equilibrium at the local temperature T(x) and mean 
velocity u^(x). The momentum distribution in each cell 
must therefore be /(x, p). The local equilibrium phase 
space distribution pO} is correlated if there are density 
correlations between cells or autocorrelations. These cor- 
relations are described by the pair density tj 2 (xi,x 2 ). 
The integral n 2 over both positions gives the number of 
pairs averaged over events (N(N — 1)). In the absence of 
correlations, n 2 (xi,x 2 ) — > ni(xi)m(x 2 ). 

In order to study the angular distribution of fluctua- 
tions, we use (18} , (fl9b, ([201), and (fl5l to write JIS) as 



KPi> P2) = / c(xi, x 2 )/(xi,pi)/(x 2 , p 2 )dTidT 2 . (21) 
We consider the angular correlation function 
Ap= r 6(A<f> -<(>! + 2 )<5(At7 - m + m)dpidp 2 , (22) 



where r(pi , p 2 ) is given by ( 10 ) . This function probes the 



(Arj, A(j>) correlations of particles in the full range of |pt|. 
Such correlations are dominated by the more abundant 
low p t particles. 

We follow Ref. [31 [33] and identify c(xi, x 2 ) with ([l6j), 
a form that describes the system at its formation. This 
identification omits the effects of diffusion described in 
Ref. |56j . This omission is reasonable only as long as 
correlations are dominated by pairs separated by |tji — 

m\ > 1. 

To clarify the contributions of fluctuations and 



anisotropic flow to correlations, we expand ( 22 ) as a 
Fourier series. Equation Q implies that the integral of 
Ap over A<fi in a rapidity range gives (N) 2 TZ. We there- 
fore write 



^ = ^{1 



2y^A w cos(nA<7J>)}. 



(23) 



Next, we compute the Fourier coefficients using (21) to- 
gether with the correlation function (15 - 17). We find 



A n = d 2 r t p FT (r t )[v n (r t )f 



where we define the local flow coefficient 



v n (r t ) = dp 



/fo.p) 

n(r t i) 



cos n(f>. 



(24) 



(25) 



Identifying the average over the probability distribution 
p FT as the event average, we se e that ( 24 ) has the form 
A n = (v%) in (KJ]); see Sec. Ill An important caveat is 



that momentum conservation corrections omitted here 
will modify the low orders, particularly A\ [44—45] . We 
will treat these corrections and present a detailed compu- 
tation of flow parameters and flow fluctuations elsewhere. 
We now compute {Spt\$Pt2)i combining (21) with (15 



17) as before to evaluate (11 ). To simplify the denomi- 
nator in (11 1, we note that (7[) implies that (N(N— 1)) = 
(N) 2 (l +K). We then obtain 



(6p tl 6p t2 ) = 



n 

l+TZ 



d 2 r t p FT (r t )[g(r t )} 



where the local momentum excess is 



g(r t ) = dp 



/(r t , P ) 
n(r t ) 



(Pt-(Pt))- 



(26) 



(27) 



In contrast to the flow coefficients (24), we see that 
{8pt\8pt 2 ) is proportional to 1Z. However, both flow and 
fluctuation quantities depend on p FT in a similar manner. 
Observe that g vanishes if the velocity and temperature 
are uniform. 



GLASMA FLUCTUATIONS 



The azimuthal dependence of the correlation function 



given by ( 23 1 , ( 24 ) , and ( 25 1 is essentially equivalent to 



that used in hydrodynamic models [19] . The relative 
height of the near-side ridge compared to the away-side 
features in A</> depends on the shape fluctuations p FT , 
the magnitude of the flow coefficients, and momentum 
conservation. However, the overall scale of correlations 
is set by the multiplicity fluctuations TZ. This scale is 
crucial if one wishes to compare correlations at different 
centralities or beam energies. The difficulty with inter- 
preting the shape of the correlation function is that one 
must disentangle information on the production mech- 
anism contained in p FT from flow and viscosity effects. 
Unless a number of simplifying assumptions hold true, 
this may prove challenging [57 . 

It is easiest to appreciate the significance of TZ in 
the context of Color Glass Condensate theory. Suppose 
that each collision produces K flux tubes, and that this 
number varies from event to event with average (K). 
Each Glasma flux tube yields an average multiplicity of 
~ aJ 1 (Q s ) gluons, where Q s is the saturation scale [55] , 
The number of gluons in a rapidity interval Ay is then 



(N) = (dN/dy)Ay ~ o.- 1 (Q.)(if). 



(28) 



In the saturation regime K is proportional to the trans- 
verse area R 2 A divided by the area per flux tube, Q~ 2 
[55] . In Ref. [32] we show that the scale of correlations is 
set by 



n 



(N 2 ) ~(N) 2 - (N) 

W 2 



oc {K)-\ 



(29) 



The dependence on (K) drops out of the product 

1ldN/dy = KaJ 1 (Q s ), (30) 

a result consistent with calculations of Dumitru et al. in 
Ref. [3TJ- Gelis, Lappi, and McLerran have shown that 
the multiplicity distribution P(N) in the Glasma follows 
a negative binomial distribution [59]. We point out that 
by definition their negative binomial parameter k NBD sat- 



isfies TZ = 1/fej 



(29) and (30) 



Their calculated k NBD agrees with 



The signature of the Glasma contribution to correla- 
tions and fluctuations is that IZdN/dy depends only on 



Q s . Equation (30) therefore constitutes a scaling rela- 
tion, since Q s depends on many collision variables in a 
combination that can be computed from first principles. 
The leading order formula in Ref. [55J relates Q s to the 
beam energy and the number of participants per unit 
area, which in turn depends on A and b. Measurements 
of the ridge at various beam energies, target masses, and 
centralities fix the dimensionless coefficient k in ( 30 ) and 



are in excellent accord with the leading-order dependence 
[32"1 [33] . Uncertainties in the underlying description of 
flow were the biggest source of uncertainty in comparing 
( 30 1 to ridge data. 



Multiplicity fluctuation measurements of TZ in princi- 
ple circumvent the complexity of flow, facilitating the 
search for this Glasma scaling or other production mech- 
As discussed in Sec. 
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FIG. 1. Prediction for IZdN/dy as a function on the number 
of participants JV P art at three beam energies (top). Calculated 
TZ compared to PHENIX data from [B0] (bottom). 



the correlation function (21), so that the cos(nA^) con- 
tributions to (23) vanish. Predictions shown in Fig. n] 



are obtained using ( 30 ) with the energy independent di- 



mensionless coefficient k fixed to fit the near side ridge 
as in (32] G3] . The number of participants N pm - t is used 
to indicate centrality. 

PHENIX has measured multiplicity fluctuations at 
RHIC [BDl HI]- They report a negative binomial pa- 
rameter k NBD . The best we can do is to compare their 

kz\„ = TZ to (301 divided by dN/dy measured elsewhere. 



Results shown m the lower panel in Fig.[T]agree fairly well 
in central collisions. 

The experimental result that TZ is not zero shows that 
geometric shape variation is not the only source of fluc- 
tuations. To deduce anything beyond that, two caveats 
are in order. First, (29) and (30) strictly apply only to 



amsm signatures. 



Ill TZ integrates 



the number of gluons. Taking the number of particles 
to be conserved through hadronization as discussed in 
Ref. [55], we can identify TZ in (29) with the measured 
multiplicity fluctuations (17]). However, one is then un- 
sure how to address phenomenological concerns such as 
resonance decay. Second, experimenters must exercise 
care in measuring TZ as a function of centrality, because 
centrality selection can distort multiplicity fluctuations. 
Using narrow multiplicity bins to select centrality will re- 
move fluctuations entirely. One way to remove this bias 
is to use a zero degree calorimeter to select centrality as 
PHENIX has done [60]. See the appendix for details. 

An alternative probe of Glasma scaling is (8pt\8pt2)- 
As with TZ, the anisotropic flow contributions vanish on 
angular integration. The quantity (Sp t iSp t 2) is designed 
to be independent of multiplicity fluctuations, reducing 
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FIG. 2. Average transverse momentum (top) and elliptic flow 
V2 (bottom) as functions of iVp ar t for the same energies as in 
Fig. n] Mean pt data in GeV are from 62-64 and V2 data are 
fromieSHBH. 
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FIG. 3. Blast-wave parameters for temperature (in GeV) and 
average velocity as functions of centrality. Data are from Ref. 

EH- 



our hadronization concerns [52 ], I53 ] , Moreover, it is ef- 
fectively free of the multiplicity bias effect, as shown in 
the appendix. Unlike 1Z, momentum fluctuations depend 
on the scale of the fluid velocity because flow enhances 
p t - We can constrain this dependence using u 2 and (p t ) 
measurements. 



In order to calculate (SptiSpt2) from (26 1 and (27 1 we 



must specify the relation between the initial transverse 
position r t of a fluid cell and its final transverse flow ve- 
locity v t . To maintain consistency with our ridge analy- 
sis in Ref. [3Sl [33J , we use a blast wave model. There we 
assumed 7V4 = Xr t and uniform temperature T, with pa- 
rameters chosen to reproduce fits to 200 GeV Au+Au p t 
spectra |32 (l33j . Sorensen et al. subsequently pointed out 
the importance of allowing for the ellipticity of the source 
in describing the centrality dependence of the ridge [H] . 
To account for the ellipticity of the collision volume, we 



now take 7V4 = \{e x xx + e y yy) 



where e 2 xy 



1 ± e. The 



average v t and freeze out temperature in 62 GeV and 
200 GeV Au+Au collisions are the same as those used 
in [321|33] and are based on an analysis in [5H]. At 2.76 
TeV the velocity is scaled up from the 200 GeV values by 
6% and the temperature is scaled up by 7%. We present 
these parameters in Fig. [3] note that the change in v t 
and T with beam energy is rather small. The eccentric- 
ity e is chosen to fit the observed centrality dependence 
of elliptic flow V2 in 62 and 200 GeV Au+Au collisions 
and 2.76 TeV Pb+Pb collisions [65H67]. The top panel 
in Fig. [2] compares these calculations to data from Refs. 
We then use ( 25 ) to compute V2 and compare 



to «2{2}, which includes fluctuations. To see how well 
these blast wave parameterizations work for the problem 
at hand, we calculate the averag e tr ansverse momentum 
(p t ) = J ptpidp/ J pidp using (18). The agreement is 
shown in Fig. [2] 

We now compute (Spt iSpg)- The multiplicity variance 
1Z is obtained from (30). We employ (26) and (271 com- 
bined with ( |17| , and use the blast wave parameters dis- 
cussed above. These are compared to data from Ref. 
[551 EH in Fig. |IJ The STAR collaboration at RHIC and 
the ALICE collaboration at LHC measure (8pti5pt2) for 
charged particles rather than pions. There is little differ- 
ence between these quantities at RHIC, but at LHC the 
K/ir and p/ir ratios are appreciably larger than expected 
in the observed range 0.15 < p t < 2 GeV [70]. Using 
the measured {p t ) and particle ratios for kaons and pro- 
tons gives (p t )ch/(pt)« ~ 1.13; PYTHIA gives ~ 1.07. 
The top solid curve in Fig. HI is our computation with the 
measured K/-K and p/ir ratios, while the dashed curve 
assumes PYTHIA ratios. Agreement with data is very 
good for central collisions where our local equilibrium as- 
sumptions are most applicable. Deviations in peripheral 
collisions may be due in part to incomplete thermaliza- 
tion, see [2"5Il2"rj] . 



VI. SUMMARY 

In this paper we have studied the connection between 
long range correlations, fluctuations, and flow. While 
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FIG. 4. Transverse momentum fluctuations {6ptidpt2)dN / dr/ 
as a function on the number of participants -/Vp art at the same 
three beam energies. Data is from [55II69] . Dashed and solid 
2.76 TeV curves represent different K and p contributions to 
charged particle fluctuations as discussed in the text. 



our results in Sec. |V| pr imarily address fluctuations, we 
discuss flow in Sec. IIVI to establish its roll in correla- 
tions so that we can isolate its effect. As anticipated, 
our flow coefficients and their fluctuations given by ( 24 ) 



and ( 25 ) depend primarily on the spatial distribution of 
flux tubes - the "shape fluctuations" - described by the 
probability distribution p FT . Consequently, the relative 
magnitude of near and away-side features in A</> correla- 
tions depends on these shape fluctuations. However, it is 
only in this relative sense that hydrodynamic response to 
initial shape fluctuations explains long range correlations 

The scale of correlations 1Z is important in comparing 
correlations at different centralities and energies. In par- 
ticular, Refs. [32J [33J compared calculations to the peak 
height of the ridge at A0 = for Au+Au and Cu+Cu col- 
lisions for a range of energies and kinematic conditions. 
Experimenters normalize the peak A/? to the number of 
pairs from mixed events ^Jp re f. When jet and momen- 
tum conservation contributions are omitted, our results 
had the form Ap/^/p^J = 1Z(dN/dy)F(A(j>), where F 
depends only on flow. We found the beam-energy, cen- 
trality, and A dependence of TZ to be consistent with 
data; see [33] for details. To relate this to the current 
context, observe that the Fourier decomposition of F 
gives ( 23 ) . This supports the arguments regarding the 



physical significance of the height of the ridge stated in 
[3TH33J S3]. A limitation of Refs. [3TH33J |43] is that it 
is difficult to distinguish flow from Glasma effects when 
concentrating exclusively on the ridge. 

Fluctuation studies provide an alternative set of ex- 
perimental techniques for attacking correlation physics. 
In Sec. |III| we showed that multiplicity and p t fluctua- 
tions are independent of anisotropic flow. Measurement 
of such fluctuations remove much of the hydrodynamic 
uncertainty of ridge studies [3TH33] [43] . We emphasize 
that the mere fact that the measured TZ and (Sp t i5p t 2) 
are nonzero proves that event- wise variation of the initial 



geometric shape is not the only source of fluctuations. 

In Sec. [V] we related the correlation strength to two 
fluctuation observables TZ and (8p t i&Pt2)- Glasma calcu- 
lations of 1Z in Fig. flland (SptiSpt2) in Fig. Hare in good 
accord with data, except for peripheral collisions. The 
deviation in peripheral collisions may reflect the break- 
down of the assumption of local equilibrium. The onset of 
thermalization in peripheral collisions modifies the corre- 
lation function. In particular, this effect has been shown 
to modify {5pt\5pt2) at low numbers of participants [25] . 
Partial thermalization describes peripheral RHIC data 
very well and, moreover, allows one to describe pp and 
AA collisions in the same model 26J. 

The fact that long range correlations can account for 
the measured (Sp t iSp t 2) leaves us with a puzzle. There 
are two types of fluctuations in nuclear collisions. First, 
each collision produces a different multi-particle sys- 
tem. Second, the partonic system produced in each 
event undergoes dynamic fluctuations as it evolves and 
hadronizes. We have focused on the source of long range 
correlations, where causality precludes the second type. 
However, fluctuation observables integrate both long and 
short range effects. Short range effects include jets and 
jet quenching, HBT, resonance decay, hadronization and 
hadro-chemical effects - not to mention novel fluctua- 
tions due to the phase transition. These effects modify 
correlations in the range A77 ^1 — 2 units or less. Jet- 
quenching, which most strongly affects central collisions, 
was estimated in Ref. [71]. Long range effects such as 
flow are charge independent. Net-charge fluctuations are 
therefore sensitive primarily to short range effects and, 
indeed, data show the appropriate rapidity dependence 
[72] . Perhaps measurement of charge-dependent fluctua- 
tions, 7Z a i, and (SptiSp t2 } a b f° r a,, b = + or — , would help 
distinguish long range from short range contributions. 
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Appendix A: Experimental Considerations 

In this appendix we discuss experimental issues that 
affect the measurement of the fluctuation observables 1Z 
and {8pti8pt2)- We show that biases can modify TZ if 
multiplicity is used to identify centrality in fluctuation 
observables, but that (Sp t iSp t 2) is largely unaffected. Ex- 
perimental considerations in fluctuation measurements 
were discussed in Ref. [30]. These observables were con- 
structed to be independent of experimental acceptance 



and efficiency effects, and many contributions that alter 
1Z were estimated. We extend the treatment in Ref. [30] 
to include (5p t iSp t 2)- 

Assume that the collision produces K independent 
'sources' which then produce particles; one can think of 
each source as a flux tube or a wounded nucleon, depend- 
ing on one's favorite model. A set of events with fixed K 
produces the single and pair densities that scale as: 

Pi=PiK, p 2 = p 2 K + p 1Pl K{K-l). (Af) 



Suppose that experimenters measure a multiplicity m 
to identify centrality. Each m bin receives contribu- 
tions from a range of K, so that the measured quan- 
tities are averaged over that ensemble. The multiplic- 
ity of a particular particle species will roughly sa tisfy 
( N )m = (J Pi) m oc (K) m . Equations 0, (|9J), and ( |A~LJ ) 
imply that multiplicity fluctuations satisfy 



n„ 



A 



(K)r 



(K 2 ) m - (K); 



(A2) 



where A is a model-dependent constant, see Ref. [30] for 
details. The first term represents the average fluctuations 
per source. The second term comes from fluctuations 
in the number of sources. Unconstrained, independent 
sources follow Poisson statistics, (K 2 ) m — (K)^ = (K) m . 
To illustrate the effect of centrality cuts on 7Z m , con- 
sider the following opposite extremes. Suppose for the 
moment that m and N are multiplicities measured in 



the same rapidity interval, so that one can approximate 
m oc K ex N. Then the fluctuations of K for fixed m 
vanish. On the other hand, if we take m to be the signal 
in a zero degree calorimeter, then m and K are correlated 
only by the impact parameter b. Fluctuations of K may 
then dominate ( A2 ) and be Poissonian [32] • 



We now consider transverse momentum fluctuations 
(8ptiSpt2)- We combine ( |A1[ ) with ( |TT| ), first noting that 
(Pt) = J PtPi/ J P\ implies that J p\5pt = 0. Therefore, 
the numerator of (lll h satisfies J p 2 SptiSpt2 oc (K) m . The 
denominator of ((ups (N (N - 1)) m = (N) 2 m (l + K m ) . 
This implies that 



(Sp tl Sp t2 ) m = 



D 



1 



(K) m i+n r , 



(A3) 



where B is another constant. This quantity has no addi- 
tional contribution from K fluctuations as in ( A2 ) . How- 



ever, a small centrality-bias effect may result from the 
lZ m in the denominator. 

To estimate this effect, observe that our CGC calcula- 
tion gives Tl ~ 0.0017 for central 2.76 GeV PbPb colli- 
sions, and 0.02 for the most peripheral value computed, 
which corresponds to A part ~ 10 participants. The ef- 
fect of centrality bias on (5p t iSp t 2) is therefore negligible. 
However, it might be important in pp collisions with a 
multiplicity trigger. Of course, the effect could be elim- 
inated by normalizing to (A) 2 rather than (N(N — 1)), 
but that would have consequences regarding the cancel- 
lation of efficiencies. 
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